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Vibrations of free and surface-coupled atomic force microscope 
cantilevers; Theory and experiment 
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With an optical interferometer, the free vibration spectra and the local vibration amplitude of four 
rectangular atomic force microscope cantilevers made of silicon have been examined 
experimentally in a spectral range of 100 kHz to 10 MHz. A good agreement with the exural wave 
theory of clastic beams was found. Coupling to torsional vibrations was also observed. When the 
sensor tip of the cantilever is in contact with a sample surface the resonances are shifted in 
frequency and the vibration amplitudes along the cantilever change. A method is presented to 
calculate this frequency shift using a linear approximation for the tip-sample interaction forces, and 
the results are compared with the frequency shift calculated from the point-mass model. The 
measured resonance frequencies of a surface-coupled cantilever do not correspond as well to the 
theoretical ones as in the free case even if the elastic-beam model is used. The reason tor the 
disagreement is found to be the geometry of the commercial cantilevers and the nonlineatity of the 
trp-sample interaction force. © 1 996 American Institute of Physics. [S0034^74$(96)02809-2] 



I. INTRODUCTION 

In the original atomic force microscope (AFM) 1 the to- 
pography of a sample surface is imaged by scanning a sharp 
sensor tip xed to a soft cantilever over the sample surface 
and by measuring the cantilever de ection. As there is a 
demand to measure surface properties such as adhesive en- 
ergy or sample stiffness with the same high Jocal resolution 
as the topography, dynamic methods such as force modula- 
tion microscopy, 2 tapping mode, 3 or atomic force acoustic 
microscopy 4 " 6 have been developed. Tn these techniques, the 
cantilever vibrates while h is scanned across the surface with 
the sensor tip in contact with the sample surface at least 
during a fraction of its vibration cycle. All these methods 
reveal changes in contrast when different materials are 
present on the sample surface. A quantitative interpretation 
of the images, however, is dif cult because two differ em 
phenomena are involved which cannot be expressed in 
simple equations: the tip-sample interaction force and the 
cantilever vibration. 

Tt is well known that an AFM cantilever is a miniature 
elastic beam with its characteristic vibration modes and 
frequencies. Resonance frequencies of higher modes of 
AFlVt cantilevers have been measured by different groups 8 " 10 
and higher vibration modes were used for imaging. 5 ' 9 ' 10 The 
local vibration amplitude of the third and fourth bending 
mode of a 23.5-mm-long rectangular cantilever made of 
glass with metal coating was measured with a beam dc ec- 
tion sensor. 1 1 The static and dynamic de ection of cantile- 
vers with triangular shape was examined experimentally and 
numerically, 3-15 and the thermal noise of a rectangular can- 
tilever was calculated considering all possible exural vibra- 
tion modes.' 6 

However, for vibration techniques which use frequencies 
around or below the rst exural vibration frequency of the 
cantilever, the elastic-beam model is commonly replaced by 
the much simpler point-mass model. As will be shown here, 
the point-mass model comes to its limits when the vibration 
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amplitude along the cantilever changes; therefore an incor- 
rect frequency shift of the fundamental mode is obtained if 
the tip-sample contact stiffness is higher than the spring 
constant of tho cantilever. A quantitative calculation of the 
resonance frequencies is important when quantitative 
information — for example, about sample surface elasticity or 
adhesion energy — is to be extracted from the images. 

II, FREE AFM CANTILEVERS 

A. Calculation of exural vibrations of free 
rectangular AFM cantilevers 

While the vihration modes of cantilevers with variable 
cross section such ag the widely used triangular cantilevers 
require numerical calculations, 12 " 15 the theory becomes sim- 
pler for cantilevers with homogeneous cross section, such as 
rectangular ones. The equation of motion for exural vibra- 
tions and its solutions for the clamped-free cantilever can be 
found in textbooks on acoustics 1 7,1 R and also in the AFM 
literature, 7 * 8 and is therefore only brie y summarized hero. 

For a homogeneous beam of uniform cross section, the 
equation of motion for exural vibrations is a differential 
equation of fourth order: 



0) 



where E is the modulus of elasticity, p is the mass density, 
A<=ab is the cross section, and / is the area moment of 
inertia, a and b are the width and the thickness of the canti- 
lever, respectively. The area moment is l=ob 3 /\2 in the 
case of a beam with rectangular cross section, x is the coor- 
dinate in the longitudinal direction of the cantilever as shown 
in Fig. 1(a), y(x) is the de ection from the rest position of 
the length element at x f fiyldx is the slope, Bf-rPy/dx 2 is 
the moment, and El- tPyldx* is the shear force at x. A gen- 
eral solution of the type 
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FJC. 1. Boundary conditions of ttn AFM cantilever with uniform cross 
section: (a) Clamped* free canrilevtr, (b) clamped spring-coupled cantilever, 
(c) clampcd-tfmicd cantilever. 



with wave number A«=2 w/\ and angular frequency 6>=2tt/ 
can be obtained by inserting Eq. (2) into Eq. (1) yielding the 
following dispersion relation: 



(3) 



a { (f- 1,2,3,4) are constants. Since the cantilever beam is of 
nhc length boundary conditions have to be fill Jlcd by 
the solutions of Eq. (3). Tf the cantilever is clamped at one 
end and free at the other end, as is the case for an AFM 
cantilever without sample contact, the de ectton and the 
slope have to vanish at the clamped end and no 

moment and shear force can be present at the free end (x 



and 



<9x 



= 0 for x = L. 



(4) 



The general solution and its derivatives arc inserted into Eq. 
(4). The resulting four equations can be solved for the coef- 
cients <?, only if the characteristic equation 



cos k a L cosh A„t + 1 -0 



(5) 



ts ful lied. The solutions {k H L, » = U,...} of this equation 
give the wave numbers k n of an in nite set of cxural vibra- 
tion modes, where n is the mode number. With the help of 
the dispersion relation Eq. (3), the resonance frequency /„ 
belonging to k n can be calculated: 
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Jn~ — ~7 — j 



with 



c c =L&? VST" 



(6) 



(7) 



c £ is a constant containing the geometrical dimensions, 
Young's modulus, and the density of the cantilever. The dc- 
ection y n (x) of each mode is described by the following 
equation; 

y»to"Jo((w» M-cosh k„x) 

cos k n L + cosh k n L \ 
~ sin kJL + swh^L (ain k "* ~ shlh k » x) } > 

where v n is the vibrational amplitude. 



(8) 



B. Comparison to the point-mass model 

In the point-mass model, the elastic beam is replaced by 
a point mass m * attached to a spring with stiffness k c . m * is 
chosen such that the resonance frequency of the system 

o>o~yfkjm* <o) 

is equal to the lowest exural vibration frequency w, . The 
spring constant k c of a rectangular cantilever is 7 

fcVa 



4L 3 



(10) 



and the lowest exural vibration frequency is obtained from 
E^s. (<S) and (7) by setting n = 1 : 



EI 

pa" 



•i**? Vt7 -0.875)=- t, 



Eb 2 
12?' 



(11) 



m* can now be calculated directly from Eqs. (10) and (11), 
yielding 



n _k c IpLab 1 



(1.875)*" 4 ^ 



(12) 



The effective mass of a rectangular cantilever is about 
one-quarter of its real mass m. This also means that the low- 
est resonance frequency of a beam with a mass m hp attached 
to its end is 



(13) 



C. Other vibration modes 

Apart from exural vibrations m the y direction [Fig. 
1 (a)l in principle, exural vibrations can also be excited in 
the direction perpendicular to .y and x (z direction). Also, 
longitudinal vibrations and Lamb modes are possible. For 
common AFM cantilevers, the resonance frequencies, how- 
ever, of theae modes are much higher man the resonance 

AFM cantilever 
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FIG. 2. Under torsional vibrations the crass section arco A of a rcctangtilm- 
cantilever with thickness b nnd width a is rotated by on angle ft 

frequencies of the exural vibrations in they direction. Ad* 
ditionally, most position detectors, such as beam dc ection 
sensors, are not sensitive to movements in the x-z plane. 
Torsional vibrations, however, cause de ections in the y di- 
rection and can have lower resonance frequencies. The equa- 
tion of motion is a differential equation of second order' 8 : 



(14) 



9 is the angle of torsion of the cross-section area (see Pig. 2\ 
C is the torsional stiffness, p is the mass density, and J is the 
polar area moment of inertia. If the cross section is rectan- 
gular, the polar moment of inertia can be calculated as fol- 
lows: 

Ton rbn t 
' Jo Jo ( * 2+J/2) dy dz= i2^ b ^ b ^^ 

(15) 

The term containing b J can he neglected because usually 
a >b. In the case a>b the torsional stiffness C is 19 



(16) 



where G is the shear modulus of the cantilever. Again, a 
harmonic solution is assumed and the boundary conditions 
for the clamped-free cantilever are inserted. Since the differ- 
entia! equation for torsional vibrations is only of second or- 
der, the phase velocity v is a constant and the modes are 
equidistant in frequency: 



2/7-1 
21 



b fd 



(17) 



(18) 



The lowest torsional vibration frequencies are about ten 
times higher than the exural frequencies for typical AFM 
cantilevers. Torsional movements of the cantilever below its 



rst resonance frequency are used in the lateral force mode, 20 
where frictional forces between tip and surface are measured. 

D. Measurement of free AFM-cantl lever vibrations 

We examined vibrations of free AFM cantilevers experi- 
mentally. This was done with the setup shown schematically 
in Fig. 3. By oscillating the clamped end, forced cantilever 
vibrations are excited. The base of the cantilever, a chip mea- 
suring 2X4X0.5 mm 3 , is xcd with coupling gel to an ultra- 
sonic transducer. The transducer ts excited with continuous 
rf in the frequency range of 100 kHz to 10 MHz. It emits 
longitudinal waves which are transmitted into the chip and 
therefore cause vibrations of the clamped cantilever end. The 
absolute amplitude of the cantilever vibrations is measured 
with an optical Michel son heterodyne interferometer. 11 The 
He-Ne laser beam of the interferometer is focused onto the 
cantilever with a microscope objective with a focal length of 
1 6.9 mm. Thus, a focal spot diameter of approximately 5 jim 
is achieved, much smaller than the cantilever width. The 
laser beam is backre ected from the cantilever and interferes 
with a reference beam in the interferometer which is 
frequency-shifted by 160 MHz relative to the signal beam. 
The interference signal is detected by a fast photodiodc, arn- 
pli ed by 60 dB and demodulated. Calibrated amplitude de- 
tection Is possible above about 100 kHz, the upper frequency 
of the phase-locked feedback loop of our demodulation 
unit To improve the signal«to-noise ratio, a low pass Iter 
of 28 MHz was used. The hered signal is displayed by a 
fast digital oscilloscope of 150 MHz bandwidth. Further- 
more, averaging of up to 200 vibration cycles and fast Fou- 
rier transforming in order to examine the frequency spectrum 
of the signals are possible. The transducer with the chip is 
mounted on a motor-driven linear translation stage with a 
smallest step size of 0.5 juun. In this way the focal spot of the 
interferometer can be scanned along the cantilever. 
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FIG. 3. Experimental setup to measure the free AFM -cantilever vibrations. 
The chip with the cantilever is ghicd to an ultrasonic transducer. The trans- 
ducer is excited with rf i n the frequency range 1 00 kHz to 1 0 MHz. It emits 
comnressional waves tvhich are coupled into the chip and csicitc forced 
cantilever vibrations by vibrations of the clamped end of the cantilever. To 
measure the cantilever Wbratiorw, the He-Nc laser beam of an optical Mich, 
elson hererodync interferometer is focused onto the cantilever. The signals 
of the inierferometer arts ampti cd and led to the input channel of a fcsi 
digital oscilloscope. With a motor driven linear translation stage, the s 
spot of the interferometer is scanned along the cantilever. 
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TABLE r. Data of the rectangular cantilevers mode of monocrystalHne sili- 
con (Rcf. 16) used in the experiments. 



Cantilever 



Length (i,) 

Width 

Thickness 

*, 

*** 

Full length (L 2 ) 



4$0 /on 
47 /mi 
2.1 jim 
0.2 N/m 
14kHz 
475 Mm 



233 Mm 
51 Mm 
1.5 Mm 
0.6 N/m 
36 kHz 
262 /Un 



230 Mm 
28 Mm 
7.4 pan 
39 N/m 
189 kHz 
245 Mm 



235 Mm 
32 Mm 
5.2 Mm 
75 N/m 
373 kflz 
142 Mm 



Four different rectangular cantilevers made of mono- 
crystalline silicon 23 were examined. The data of the cantile- 
vers are given in Tabic I. Length L , , width, thickness, spring 
constant, and resonance frequency were measured or calcu- 
lated by the manufacturer. L } is the length from the mount to 
the sensor tip. As can he seen from the optical micrographs 
(Fig. 4), the full length L 2 of the cantilevers is slightly larger. 
Cantilevers 1 and 2 are relatively soft and thin compared to 
cantilevers 3 and 4. In all cantilevers the (110) direction is 
parallel to the length axis. The Young's modulus E of silicon 
in this direction is 1.69X10 1 ' N/m 2 (Ref, 24) and the mass 
density p is 2330 kg/m 3 . 

In a rst experiment, the amplitude of the vibration 
modes was measured along the cantilever and compared with 
the theoretical values given by Eq. (5). Jn Figs. 5(a)-5(c) the 
squares and the solid line show the vibration amplitude mea- 
sured with the interferometer for cantilever 2, The dashed 



line represents the calculated amplitude. 1=260 fjun was 
assumed. In every graph, the amplitude was ttcd using only 
one point on the graph. The clamped end of the cantilever is 
at jc = 0\ The mode /r= 1 was not measured because its fre- 
quency is smaller than 100 kHz and hence below the cutoff 
frequency of the interferometer. A very good agreement be- 
tween theory and experiment can be observed. Near the free 
cantilever end, the signat-to-noisc ratio of the measured sig- 
nal decreased because the laser spot was refracted by the 
sensor tip or by the edges of the cantilever. 

From our measurements it becomes clear that the vibra- 
tion Spectrum of a cantilever cannot be measured without 
scanning the detection spot along the cantilever. Otherwise 
one might miss a mode because h accidentally might have a 
node at the spot position. The only position on the cantilever 
where all modes have a maximum of their vibration ampli- 
tude is the free end. 

The exuraJ vibration spectra of the four cantilevers arc 
shown tn Figs. 6(a) -6(d). The amplitude was measured in 
the maximum nearest to the free end The excitation fre- 
quency was changed stepwise manually, and afterwards the 
spot position was readjusted slightly to nd the maximum of 
the vibration amplitude. The modes with frequencies below 
100 kHz could not be measured as mentioned above. For all 
cantilevers the spacing of the modes increases with increas- 
ing frequency and with increasing cantilever stiffness in ac- 
cordance with theory. In Table II the experimental exural 
resonance frequencies/^ are compared to the calculated f n . 
The corresponding wavelength \ n of the modes is also given. 




I 1 

100 pm 




(b) (d) 

FIG. 4. Optical mrcrt^pha shmvlng (a) cantilever 1. (b) cantilever X (c) contikra 3. and (d) cantilever 4. 
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(a)n = 2, 

f = 191.9 kHz 




(b) n = 3, 

I = 547.3 kHz 




x [tun] 



(c) n * 4, 

f = 1.082 Mhz 




<d) 

f ^ 1.705 MHz 




(e) n b 6, 

f =2,sea MHz 




x [Mm] 



FTO. 5. Vibration amplitude or mode /i=2-6 along cantilever 2 measured 
wfth the optical interferometer. The clomped end is at jt-0 and the free end 
at x - 1. The dashed line shows the theoretical amplitude v„(x) ttcd to the 
experimental values using one point on every graph. A length or 260 fun 
was assumed for the calculation. 



For all cxiiraJ modes 2 should give the constant 

value Uc\ according to Eqs. (6) and (7). 

Figure 7 shows the difference of /^(^L) 2 from the 
average value in % for the four cantilevers. For cantilever 1 
/oii/(*ji^) 2 is almost perfectly a constant; only mode n= 12 
obviously does not t into the cxural vibration spectrum. 
The vibrational spectrum of cantilever 2 also shows a good 
agreement between experiment and theory. The mode /i = 8 
could not be excited experimentally. The two stiff cantilevers 
3 and 4, however, show an appreciable deviation from 
theory. 

The mass of the sensor tip was not considered in the 
calculations- The sensor tip is a small cone of J 5 fjem height 
and 35° opening angle. Its mass is approximately 
8X10" 13 kg. This means that for all cantilevers the effective 
mass m*«-/7j/4 is more than one order of magnitude larger 
than m^. 

Since in the elementary theory of cxural vibration ro- 
tation and shear de ection of the cross sections arc not taken 
into account, the phase velocity increases in nitely with in- 
creasing mode number. A theoretical description taking both 
effects into account leads to a phase velocity converging to- 
wards a value between the Rayleigh wave velocity and the 
bulk shear wave velocity. 25 However, as long as the cantile- 
ver thickness divided by the mode wavelength b/\<0A 9 the 
elementary theory is in good agreement with the exact 
theory- From Tables .1 and 11 one can see that this holds for 
all cantilever modes, even for the stiff cantilevers. The sys- 
tematic deviation of the resonance frequencies of the stiff 
cantilevers is due to the cantilever shape. At the tree end, no 



uniform cross section is realized as assumed. Also, the can- 
tilevers arc not real ly clamped at x = 0 but continued as a sort 
of stripe waveguide on top of the chip which changes the 
boundary condition. 

From the width of the resonance peaks A/, the quality 
factor e=/ wp /A/ was calculated for all resonances and is 
shown in Fig. 8. For the cantilevers 1-3 Q has a maximum 
between ) and 2 MHz. The damping of free AFM cantilevers 
is mainly caused by sound radiation and friction with air. In 
vacuum, the Q values of the AFM cantilevers are orders of 
magnitude higher than in air. 8 

As already discussed earlier, torsional vibrations must 
also be taken into account when analyzing AFM cantilever 
vibrations. If the cantilevers are slightly asymmetric, cou- 
pling between torsional and exural vibrations becomes pos- 
sible. Mode coupling between a cxural and a torsional 
mode probably happened when we tried to excite mode n = 8 
of cantilever 2. When the excitation frequency was tuned to 
5.05 MHz, the cantilever started to vibrate at the two neigh- 
boring frequencies of 2.35 and 2.7 MHz, whose sum is 5.06 
MHz, instead of vibrating at its resonance around 5 MHz. 
The peak at 2.7 MHz can be assigned to the exurai mode 
« = 6 which coupled to a torsional mode of 2.36 MHz reso- 
nance frequency. Also, the vibrational mode with the twelfth 
frequency which docs not t into the spectrum of cantilever 
1 is probably a torsional mode. Since the torsional modes arc 
equidistant in frequency while the spacing of the cxural 
modes increases, the probability increases to excite a tor- 
sional mode at higher frequencies. 
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FTG. 6jTcqucncy S pcctm of (a) eanjilcvcr 1, <b) cantilever 2, (c) cantilever 3, mid (d) cantilever 4 measured experimentally with the optica] interferometer. 

J^l^T^l^ ^ , m thC TCR0 " 3nce P^ 9 " Mo « ° f * c Tcsonance peaks appear » vertical lines because their nilc width is not resolved 
in ine ircqucncy scale used to display oil measured resonances. 



The comparison between experiment and theory for the 
free cantilever has shown that the exural vibration theory is 
suf cient to describe AFM cantilever vibrations in the fre- 
quency range we discussed here. Nevertheless, the geometry 
of the cantilevers plays an important role, and for the stiffer 
cantilevers a better agreement with experiment can probably 
only be achieved by calculating the modes numerically. In 
the special case regarded here, the mass of the tip could be 
neglected because it is very small compared to the cantilever 
mass. The rst four to ten modes could be excited with high 
amplitude, with Q factors ranging from 90 to 905. This 
means that the higher modes arc well suited for sensing pur- 
poses. Using the higher modes with many nodes, special care 



has to be taken in order to choose the position along the 
cantilever on which the vibration amplitude is measured. 

III. AFM CANTILEVER WITH ITS SENSOR TIP 
CONTACTING A SAMPLE SURFACE 

A. Linear model for the cantilever vibrating with 
surface contact 

When the cantilever is approached to a sample surface, 
in most cases the sensor tip senses rst the long-range attrac- 
tive tip-sample interaction forces such as Van der Waals 
forces. Using a soft cantilever, the tip jumps into contact 
with the sample surface into an equilibrium position, where 
the tip-sample interaction force is balanced by the cantilever 



Cantilever I 



Cantilever 1 



Cantilever 3 



(MHz) 



(MHz) 



Cantilever 4 



(MTU) 



(MHst) 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 



1.875 
4.694 
7.R85 
10.996 
14.137 
17.279 
20.420 
23.562 
26.704 
29.845 
32.987 
36.128 



0.22383 

0.4395 

0.7279 
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2.021 
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3.94 

4.24 



0.013 
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1.0S4 
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2.016 
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4.741 



1592 
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146 
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3 00 
90 
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0.031 
0.195 
0.550 
1.069 
1.767 
2.640 
3.686 
4.908 
6,305 



878 
351 
209 
150 
116 
95 
81 
70 
62 



(MHz) 


f* 


X„ 


few 

(MHz) 


(MHz) 


(A*n>) 


0,16305 


0.157 


831 


0J3085 


0.310 


476 


1.0052 


0.984 


332 


K986 


1.943 


190 


2.762 


2.775 


198 


5.337 


5.4R2 


113 


S.265 


5.398 


142 


9.99 


10.662 


si 


8.64 


8.922 


no 







3288 



Rev. Scl. Instrum., Vol. G7, No. 9. September 1096 ACM m 

Aral cantilever 

PAGE * RCVD AT 3/1 1/2005 7:40:09 PM [Eastem Standard Time] * SVR:USPT0^XRF-1/1 * DmS:87M306 * CSID:4083820481 * DU1RATI0H (mm-ss):1 1 ^0 



03/11/2085 .16:49 4083820481 



PATENT LAW GROUP LLP 



PAGE 21/27 




H H f- 

9 10 11 12 



AG. 7. Difference A in % or/^A/) 2 fVom its average value ibr the Pour 
cantilevers examined experf mentally. For cxural vibration modes 
fA k J<) 2 should he a constant. This is the case fhr most of the resonance 
frequencies found for cantilevers I and 2 while the stifTbr cantilevers 3 and 
4 show a deviation from this behavior. 

force. When an ultrasonic wave is coupled into the sample, 
such that the sample surface vibrates at ultrasonic frequen- 
cies, the sensor tip and the cantilever are also forced to vi- 
brate around their equilibrium position. The tip-sample in- 
teraction forces are highly nonlinear. However, if only very 
small vibration amplitudes of the cantilever around its equi- 
librium position are excited, the tip-^amplc force can be ex- 
panded linearly, and H can be represented by a linear spring 
with spring constant **, where k* is the negative derivative 
of the tip-sample fbtec in the equilibrium position: 26 



09) 



Here f is the tip-sample distance, is the tip-sample 
interaction force, and | c is the equilibrium position. 

The spring model is shown in Fig. 1(b). In the contact 
case, the forces exerted by the spring /c* arc repulsive when 
the cantilever approaches the surface and attractive when it 
moves away from the surface. As can be seen in Fig. 1(b), it 



is assumed that the AFTvl cantilever is positioned parallel to 
the sample surface, as a consequence the exural cantilever 
dc cctions are normal to the sample surface., i.e., in the £ 
direction. Furthermore, it is assumed that the forces which 
act onto the sensor tip act as a point force on the free end of 
the cantilever. In the simple model used here, the cantilever 
is not de ected in its rest position. In the contact mode, the 
corresponding equilibrium position of the sensor tip can al- 
ways he adjusted experimentally by moving the sample sur- 
face relative to the cantilever mount. 

To analyze the behavior of this system, the resonance 
frequencies of a cantilever with one clamped end and the 
other end coupled to a spring must be calculated. The spring 
at the end changes the boundary conditions at x=*L. Like the 
free end the spring-coupled end can transfer no moments. 
But a de ection y generates a force -k*y which has to be 
added to the shear force EltiPyldx*: 

*/^T-A*v=0. (20) 
Therefore, the boundary conditions at jr m L are now 



*y 
J? 



= 0, 



3A* 



\ y for x=L. 



(21) 



The boundary conditions at x= 0 are the same as in Eq. (4). 
Instead of using Eq. (2) it is more convenient for the calcu- 
lation to express the general solution for the local amplitude 
(omitting the time dependence) in the following way; 

v(j:)=>J(cos Ajc + cosh *jc)+£(cos Ax-cosh kx) 

+ C(sin to+ainh fcc)+£>(sm Jbc-sinh kx). 

(22) 

The characteristic equation for this system is calculated in 
the same way as for the clamped-free cantilever by inserting 
the general solution and its derivatives into the boundary 
conditions and solving for the constants A, 5, C„ and D in 
the general solution. After some straightforward calculation 
one obtains 
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TABLE 1IL Boundary conditio** al x-L and characteristic equations tor the throe different CoScs shown in 
Fjjs(& 1 Co)- 1(c). 



Boundory conditions 
at * = 



Characterl?lk equation 



(o) Free end 
(in nitely soft 
spring. A*~0) 




*> • 


cos a^x, co.tn Kjflt + 1=0 


(b) Spring, 
coupled end 
(contact mode, 
repulsive tip- 
sample forces) 


a? 


* * 
*• 3** 


Jrinh k„L cos k n L -sin cosh k H L 


(c) Pinned end 
(in nitely saJfT 
spring. k m =*>) 




5?"° 


sinhA,/. eos* fl L— sin *„I cash * n £»0 



sinh k„L cos A„Z. - 9in cosh k„L 
= ~ 3p — (1 + cos cosh 



(23) 



Equation (23) reduces to Eq. (5) if ** = 0\ i.e., if the 
spring is in nitely soft In the other extreme case, when 
= one obtains a cantilever pinned at one end [Fig. 
1(c)]. The boundary conditions at * = Z, and the character- 
istic equation for the damped-pinncd cantilever arc listed in 
Table III. The solutions of Eq. (23) were found numerically 
using Mathematica} 1 Some examples are listed in Tabic JV. 
In Pig. 9 the rst three resonance frequencies f n of a spring- 
coupled cantilever are shown as a function of A*/A c ./„ was 
normalized to the rst resonance frequency of the clamped- 
free cantilever. The resonance frequency of every mode of 
the clamped-spring-coupled cantilever is between the reso- 
nance frequency of the clamped-frec cantilever mode and the 
resonance frequency of the clamped-pfnned cantilever mode. 
The open circles in Fig. 9 show the resonance frequency f n% 
of the spring-coupJed system normalized to the free reso- 
nance frequency f Q as it is calculated from the point-mass 
model: 36 



fo " 



4 



(24) 



If the tip- sample contact is softer than the cantilever 
spring constant, fc*<k Cf there is a good agreement between 
the point-mass model and the elastic beam model. However, 
if k*>k €t k c becomes negligible compared to k* in the 
point-mass mode). This means that the system now behaves 
like a poin t mass m* attached to a spring k* and f m 
>lk*/m*. The elastic beam, however, can vibrate even if 
both ends are xcd. Therefore the point-mass model fails if 
k*>k c . The local vibration amplitudes of modes 1, 2, and 3 
can be seen in Fig, 10. They were obtained using Eq. (8) and 
the k n L calculated from Eq. (23). For very soft springs the 
local vibration amplitudes are maximal at the end of the 
cantilever and increase compared to the clamped-frce case. 
When the spring becomes sttffer, the slopes at the end of the 
cantilever go through zero and change their sign. If A* in- 
creases further, the maximum of the vibration amplitudes is 
no longer at the end of the cantilever and the vibration am- 
plitudes at the end decrease until they become zero in the 
clamped-pinned case. Note that the k* which causes zero 
slope at the end increases with increasing mode number, i.e., 



Free end 
(*•=[» 



= 0.1 



1 
2 
3 
4 
5 
6 
7 

a 



1.8751 
4.69409 
7.85476 
10.9955 
14.1372 
17.2788 
20.420$ 
23.5619 



1.91891 
4.69699 
7.85538 
10.9958 
14.1373 



2.2135 
4.7234 
7.86097 
10.9978 
14.1382 
17.2793 



3.16765 
5.00)12 
7.91896 
11.0185 
14.1479 
17.2846 
20.4239 





**/*,= I0D0 


Pinned end 
(**-«) 


3.82981 


3.91704 


3.9266 


6.40415 


7.00766 


7.0685B 


8.58742 


10.012 


10.2102 


11.2559 


(2.851 


13J518 


14.2523 


15.4765 


16.4934 


17.3396 


18.0795 


19.635 


20.4567 


2DJT774 


217765 


23.5853 


23.839 


25.9181 
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FIG. 9. Resonance frequencies /„ or Ihc clamped spring<oupled system 
(black squares) normalized to the rst resonance frequency Of the clamped- 
free cantilever / 0 . As shown here for Ihc rst three modes, ihc resonance 
frequency of ihc nth mode of the clamped npring-eouplcd cantilever lies 
between the frequency of the wth clamped- free mode and the nth damped- 
pinned mode. When **/* c « increased from 0 to * the resonance rrccrucney 
of every mode shifts from the free case to the clumped case. The gray area* 
indicate frequency gaps. A comparison with the point-mass model for the 
cantilever (open circles) shows that this model predicts too large frequcncY 
Shite for A*/A C >1. 



a very soft spring has little in uencc on the higher modes, 
while a very stiff spring cau9cs the lower modes to behave 
almost like in the pinned case. 

Additionally to a frequency shift of the resonances one 
observes an increase in the damping when the sensor tip 
comes close to the sample surface. There also might be a 



Mass m H 




- -T Dash pot y 

^^^^ 



FIG. 1 1. Additionally to the forces between the sample surface and the 
sensor tip— represented by the spring, damping-^ represented by the 
dashnot- -has to be considered. Also there mlRhi be a mass tu^ xed to Ihc 
end of the Cantilever. 



mass attached to the end of the cantilever. Both effects 
result in a modi cation of the forces at x~L (Fig. 11) and 
therefore change the fourth boundary condition: 



Ac' 



=0, 



dp 



(25) 

y is the damping constant, and a damping force proportional 
to the relative tip-sample velocity was assumed. The time 
derivatives in Eq. (25) can be calculated easily, and therefore 
the last boundary condition can be expressed as 



Ac 3 "" 17 y( k * + 'rw-WtipW 2 ) 



(26) 



or, using Eq. (3), 



n = 1 



n = Z 



n = 3 




0 0.2 0.4 0.6 0.8 1 
X/L 



0 0.2 0.4 0.6 0.8 1 
x/L 



0 0.2 0.4 0.8 0.8 1 
x/L 



undo eet«rf rest potions when the end of tte ^c^fcT^i M In * LIZ^^T Z "'"V" <m top. Tlx springs an nlways fa ft* 
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a?" 



with 



(27) 



The characteristic equation is then similar to Eq. (23): 
Rtnh k„L cos k„L-sin k n L cosh kj. 
kl 



<p(k) 



(l + co<* k n L cosh k„L). 



(28) 

This equation also gives the exact solutions for the free vi- 
bration frequencies (** = 0, y=0) of a cantilever with a mass 
at its end. In this case <p(k) reduces to 
<PW - - k^m^pA = - k A Lm bp /m and the characteristic 
equation is 

1 +cos k„L cosh k n L = - . A B L(smh A n £ cos *„Z 

-sin *„L cosh *„/:). (2.9) 

To solve Eq. (29) only for the lowest wave number of 
the system, i.e., for small ft, the trtgnonornetric functions can 
be expanded in polynomial series: 

iMD 2 (M^ 

41 



1+1 



21 + 4! j\ A ' 2! 



1 



31 



(30) 



When all fectors containing k x L higher than the fourth power 
are neglected, the approximation for the lowest resonance 
frequency reduces to 



m 



(31) 



This approximation is equivalent to Eq. (13) as can be seen 
when (1.875)*M2 is inserted and k x L is calculated from 

B. Cantilever vibrating in contact with a sample 
surface experiment 

When the sensor tip contacts the sample surface, forced 
vibrations can be excited by sample surface vibrations. To 
examine the spectra of surface-coupled cantilevers experi- 
mcntaHv we used our atomic force acoustic microscope 
setup. • A piezoelectric transducer was bonded to the 
sample with a eoupling gel. We excited the transducer with 
an electrical spilce, so that the sample surface movement was 
a short pulse containing a broad band of frequencies. The 
resulting cantilever vibration was measured with a knife- 
edge detector 29 and Fourier transformed. The sample was a 
thin glass plate partially covered with chromium. At high 

32*0 Rev. Sd. mstrunu. vol. 67. No. ft, September 1996 
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HO. 12, Frequency spectra of cantilever 3 measured on top or u gloss 
saroprc partially covered with chromium. The dashed lines show the i 
nance frequencies when (he cantilever is free. 



excitation voltages the cantilever vibration spectra contained 
many signals, but at very low exciting voltage most of the 
peaks disappeared and the remaining resonances could be 
assigned to the theoretical modes. 

Figure 12 shows the spectra which were obtained with 
cantilever 3. The dashed lines are the free resonance frequen- 
cies which were measured in the experiment described in the 
previous section. In accordance with theory, the frequencies 
of all modes are higher in the contact case. In Table V(a) the 
resonance frequencies arc listed The frequency shift is larger 
on the chromium-covered parts of the sample than on the 
glass surface. Using Eqs. (6), (7), and (23), in principle A* 
can now be calculated from the shifted resonance frequencies 
in the spectra: 

TABLE V. Resonance frequencies of the spectra shown in Fig. 12 qnd 
corresponding **. 

(a) Experimental resonance frequency or cantilever 3 on top or the glass 
sample partially covered with chromium. An accuracy of ±5 kHz is as- 
sumed 



Mode n 



I 



Frequency glass 
Frequency chromium 



675 kHz 
725 kHz 



1.825 MHz 
2.075 MHz 



3.175 MHz 
3.675 MHz 



(b) k+ calculated with the data given by the manufacturer: c r ~ 4.274 
X lO^vfe, modeled with spring al the end. The errors were calculated by 
varying the resonance frequencies by ±5 kHz 

Modon ] 2 3 



k* (N/m) gl as5 
k* (N/m) chromium 



863+30 
I293±60 



1617±17 
2703 ±28 



-1850 ±60 
1863±29 



(c) k* calculated with the data derived from the free vibration frequencies. 
Moden = l:c r * (4.6437:* 0.001) X 10"' & „=2: e B - (4.682 ± 0.0O1) 
x I0"-' > /s:« = 3;c < . = (4.7263 ± 0.0005) X 1 0 " 3 v5 

Mode n ] 2 3 

** (N/rn) glass 2644 ±350 2621139 2234 ±2B 

ft* (N/m) chromium - 1 0527±8000 fi046± 1 65 5 1 75 ±4) 

(d) A # caIculotcdwftht,=(230±5)^m(c, = (4.274 ± 0.1) X Ur*f*), 
and i 2 derived from the free resonance frequencies as in (c) 

Mode n 1 2 3 



*• (N/m) 



1005+130 
1428*370 



3345±410 
471 3 ±440 



9198±3870 
13901 ±5040 
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x'=0 



length L ' extending from the spring to the free end. The full 
length of the cantilever is L^Li+L'. The characteristic 
equation of this System is 

-*(cosh k n L x sin k n L { ~s\nh k n L^ cos k n L{) 

x{l + cos k rt L' cosh * n £') + (cosh k n L' sin k n V 

-sinh k n L* cos * n £')(l-ccs k„L\ cosh kj, x ) 

^^.^^[J+cos^U^L^coshAjA^+L')]. (33) 

A derivation of this equation is given in the Appendix. The 
results for ** obtained with this equation using £, = 230:3:5 
Mm and the free resonance frequencies to determine L 2 arc 
shown in Table V(c). Now, no resonance frequency is in a 
gap, however, k* increases with frequency for both, glass 
and chromium. 



FIG. 13. Mode) fhr condJevw* with a sensor tip xcd at any lentfh dement 
of the beam. The cantilever can be divided into two parts with length I , and 
L' and individual solution* ^(jc) and /(*'), respectively. As in the simpler 
case with ihc spring at the end as shown in FJga. |(a)«-l(c), the spring- 
coupled boundary condition converges to the pinned boundary condition 
when **—»«=. 



k* = k 



l+cos(c cV ^)co R h(c cV ^) 



e Sinh(c c >//;)cos(c c y[f n )- stn(c c v^)cosh(c c ^) 



( 32 > 

Table V(b) shows the results which were obtained for k* 
using the cantilever dimensions and k c given by the manu- 
facturer and in Table V(c) A* was calculated using the free 
resonance frequencies measured experimentally to determine 
c c . This comparison shows that the values obtained for k* 
depend very sensitively on the geometrical dimensions of the 
cantilever inserted into Eqs. (7) and (23). As can be seen in 
Fig. 9, there arc gaps in the frequency spectrum of an indi* 
vidua! cantilever, t.e., frequencies which cannot be resonance 
frequencies of the clamped/spring-coupied system. Negative 
values for k* are obtained when the experimental resonance 
frequency is in such a gap. The large errors in Table V( c ) for 
mode n=\ on chromium arc obtained because the tip- 
sample contact stiffness is much higher than the spring con- 
stant of the cantilever. Therefore, the frequencies of the low- 
est cxural vibration mode already converge towards the 
frequency of the clamped-pinned case, which means that 
Amal] changes in resonance frequency correspond to very 
large changes in k*. As pointed out above, a disagreement 
between the cxural wave theory and the experiment could 
be observed originating from the inhornogencous cross sec- 
tion of the cantilever and from the boundary conditions at the 
clamped end. The fact that the sensor tip is not exactly at the 
end of the beam (sec Fig. 4) can also cause frequency shift? 
which do not agree with the simple theory and a situation 
arises as shown in Fig. 13(a). Now the cantilever can be 
divided into two parts, the rst part with length L } extending 
from the clamped end to the spring and the second part with 



C. Calculation of K* from Hertzian contact theory 

To compare theory and experiment also for surface- 
coupled vibrations, the tip-sample interaction force must be 
known in order to calculate A*. Especially when working in 
air, these forces are very complex and cannot be described in 
a simple way. A rough estimation of the order of magnitude 
for** can be obtained by using the formula for the Hertzian 
contact theory 30 when the tip penetrates into the sample sur- 
face, In the experiment the force applied by the cantilever 
was adjusted to be negligible compared to the adhesion 
force. AH attractive adhesion forces arc summed in a force 
F Q which is treated like an external force applied onto the 
tip. F 0 can be estimated from the so-called force calibration 
curves, where the tip is lowered until it is attracted by the 
sample surface and afterwards withdrawn until it detaches. 
The Hertzian contact force is 



/ T- 



with 



_1 \-v 2 ( l-v] 



(34) 



(35) 



where R is the radius of the sensor tip and E* is the effective 
Young's modulus of the tip-sample contact. E t , v n E & , and 
v a arc the Young* s moduli and the Poisson ratio of the tip 
and the surface, respectively. The equilibrium position g e is 
calculated from 



4£* i 1— (3F«\™ 



(36) 



Calculating the derivative of Eq. (34) and inserting Eq. (36) 
yields 



= -\j6E* 2 -R.F 0 



't-tr. 

and this means that 
A»=V6 E^.R-Fo. 



(37) 



(38) 
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The sensor tip is oriented in the (100) direction. The 
Young's modulus and Poisson ratio of silicon in this direc- 
tion are 1.3098xJO n N/m 2 and 0J.81, respectively. 31 The 
Young's modulus of glass range* from (0.4-0.9)X 10 11 N/m* 
and the Poisson ratio is O.2.* 2 For the Young** modulus of 
chromium such different values as 1.32 and 2.53X 10 1 1 N/m 2 
can be found in the literature. 32,33 The Poisson ratio is 0.3. 
The effective elasticity E* for the sensor tip on the glass 
surface is therefore E* = (3.2-5.5) X 10 10 N/m 2 and the ef- 
fective elasticity for the sensor tip on chromium is 
(7.0-9.1)XIO JO N/m 2 . To calculate k* the sensor tip radius 
H and the adhesion forced are needed. The radius of cur- 
vature of fresh sensor tips used in our experiments is 10 nm. 
However, the radius probably became larger when the tips 
were used for some time, especially when they were vibrat- 
ing on a stiff sample surface. From the smallest features 
which can be imaged with our tips, however, it can be con- 
cluded that R cannot be larger than 100 nm. F 0 varies in our 
experiments in air typically from 5xl0" 8 to 5X10" 7 N. 
When all the uncertainties are considered, k* can range from 
145 to 973 N/m on glass and from 245 to 1355 N/m on 
chromium according to Hertzian contact theory. 

Comparing these results from Hertzian contact thoory to 
the experimental stiffness in Tables V(b)-V(d), it can be 
seen that the result for mode »=l in Table V(b) agrees 
roughly with the upper limit for k* calculated from Hertzian 
theory. However, the other two modes yield different results, 
even one with a negative value. Using the constants derived 
from the free resonance frequencies, one obtains roughly the 
same k* for all three modes, as shown in Table V(c). Only 
on chromium the resonance frequency of mode n — 1 yields a 
negative value for k*. For glass and chromium the experi- 
mental k* in Table V(c) is more than two times higher than 
the upper limit from Hertzian theory. If one takes into ac- 
count that the sensor tip is not at the end of the cantilever. A* 
again becomes smaller and therefore closer to the theoretical 
values at /> = ! [Table V(d)], but now k* increases with 
mode number (or frequency). 

IV. DISCUSSION 

The comparison has shown that the frequency shift of 
the resonances caused by the tip-sample contact elasticity 
depends very much on the geometry of the cantilever. Be- 
cause the geometry of the cantilever plays such an important 
role for the determination of elastic constants of a sample 
surface, cantilevers with a simpler and better dc ned geom- 
etry are needed. It is important that a stiff clamp is realized 
also for the stiffer cantilevers, that the beam cross section is 
uniform, and that the sensor tip is at the end of the cantilever. 
A numerical characterization would be too time consuming 
because of the limited life span of each cantilever. 

On the other hand, it is very possible that the stiffness of 
the upper surface, layer penetrated by the sensor tip is higher 
than the values calculated from the bulk elastic moduli of the 
materials. In further experiments and calculations, the in u- 
ence of damping forces will be examined. More precise con- 
tact mechanical models considering adhesion forces have to 
be used. However, it is possible that macroscopic models 
which treat the sensor tip as an homogeneous elastic body 
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and which consider no plastic deformation and material 
transfer between tip and surface do not describe the contact 
of a sensor tip of a few nanometers of radius accurately. One 
main problem in atomic force microscopy is the unknown 
shape of the sensor tip. From changes in image contrast we 
conclude that sensor tips change their shape when they are 
used in the contact mode. Stiff materials such as silicon are 
more likely to break than to be plastically deformed. If, for 
example, the silicon sensor tip used in the ultrasonic contact 
mode was broken, the spherical shape of the sensor tip was 
lost and the high measured contact stiffness might have been 
caused by a larger contact area than the one which was cal- 
culated from the assumed spherical shape of the tip. 

Furthermore it has to be kept tn mind that knife-edge and 
bcam-de ectton detectors measure the local de ection and 
not the amplitude. A foil characterization of the cantilever 
vibration in contact with a sample surface could only be 
obtained by scanning the cantilever with an interferometer 
similar to the experiments we performed with free cantile- 
vers in order to measure the calibrated local vibration ampli- 
tude along the x axis. The accuracy in the measurcrncrrt of 
the contact resonances can be improved by exciting narrow 
band signals at well de ned frequencies instead of Fourier 
transforming only one pulse containing a broad band of fre- 
quencies. 
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APPENDIX: CALCULATION OF THE CLAMPED/ 
SPRING-COUPLED SYSTEM WITH THE SPRING FIXED 
AT A GIVEN LOCATION ALONG THE CANTILEVER 

As shown in Figs. 13(a) and 13(b) the cantilever can be 
divided into two parts. The ret part extends from the 
clamped end at x=* 0 to the spring at x =£, ; the second part 
extends from the free end at x' = 0 to the spring at x = L\ 
The full length of the cantilever is L 2 ^L l +L'. The origin of 
the x* axis is chosen to be at the free end to simplify the 
calculation. The vibration amplitude is labeled y(x) for the 
rst part and y'(x l ) for the second part, respectively. The 
boundary conditions for^jt) my{0) = 0 and <?y/dx{0) = 0 
at the clamped end and the boundary conditions for y'{x') 
are *y/<fc'*(0) = 0 and ^y7**'*(0)«0 al the free end. 
Solutions for the whole cantilever can be found from tting 
y(x) and v'(jf') together at the length element x=£, and 
hence x' = L \ taking into account the boundary conditions. 
To calculate the characteristic equation in this way is 
straightforward but a tedious task and therefore omitted here. 
However, it is much simpler to calculate the characteristic 
equation for the situation shown in Fig. 13(b), where the 
spring is replaced by a pin, and where y{L x ) = y' (£') = 0 at 
the pin. In this case, 

j>(*)=*(cos A*-cosh **)+Z>(sin Ax-sinh kx), 

(Al) 

AFM cantilever 
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jr'(x'W(cot Ax' + cosh kx') 

+ C"(sin fo'+sinh kx') (A2) 

with 

B sin ALj — sinh kL^ 

0 cos | - coah A£ i 

and 

A' __ 3in*L' + sinh kV 
C ~ cos^'-f cosh kL" 

At the length element where x = L } and;e' = £' the following 
boundary conditions have to be ful lied; 

= U,> — £(!/>■ (A3) 

5?Ui)-i^(t'). (A4) 

The third derivatives need not be equal because the dif- 
ference of the forces exerted by the two pans of the cantile- 
ver beam arc carried by the pin. When the derivatives of Eqs. 
(Al) and (A2) are inserted into Et|S. (A3) and (A4), one 
obtains the characteristic equation of a cantilever pinned at 
any position along its length axis: 

(cosh h n L\ sin * ff £,-sinh k rt L } cos *„Ll)0 + cos k„L f 

Xcosh (coah k„L $ sin k a L' 

-sinh*„Z' cojA^'Xl-cos*^! cosh k n L } ) = 0. 

(AS) 

In Eq. (23), the characteristic equation of the clamped- 
pmned case and the characteristic equation of the clamped- 
free case arc coupled by the factor = ( A„£ ) 3 * r /3&*. 

In the situation treated now when A*— >0, the cantilever vi- 
brates like a clamped-free cantilever with length L 2 = T, K +Z/ 
and the wave numbers k„ of the cantilever are given by 

1 + cos kJLt + */)cosh k„lL } +Z ') - 0. (A6) 

The characteristic equation for the clamped-pinncd case ts 
given in Eq. (A5). This means that the characteristic equation 
for the clamped-spring-couplcd case is analogous to Ea. 
(23): 

(cosh k n L { sin k n L^9\nh k n L { cos k n L x ) 
X(l+cos*„L' cosh kj,') 
-(cosh k n L 9 sin k n L' -smh k„L' cos k n L>) 
X(l -qos k n L } cosh k n L{) 
% EI 

=A-k\ ^ [1 +cos *„(£, +//)cosh k n {L x +£')]. 

(AT) 



The value of the factor A can be found by considering the 
case l'— >0, i.e., if the spring is moved to the end of the 
cantilever. Only if /f=-2, Eq. (A7) is equal to Eq. (23). 
When damping and the mass of the sensor tip are to be 
considered, k*!EJ can be replaced by <p(k) as in Eqs. (27) 
and (28). 
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